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Abstract 

Tile  anaiysij  of  eye-movemeni  (EM)  signet  pose*  problems 
for  (Jie  designer  of  smoothing  filters  since  many  of  the  in¬ 
teresting  types  of  EMs  are  bimodsJ.  For  example,  optoki¬ 
netic  and/or  vestibular  stimulation  results  in  an  EM  pat¬ 
tern  called  nystagmus  consisting  of  alternating  fast-  and 
slow-phase  components.  Also,  saccadic  (refixation)  EMs  do 
not  occur  continuously,  but  are  interspersed  with  periods 
of  fixation.  Conventional  linear,  low-pass  filters  (both  fi¬ 
nite  impulse  response  (FIR)  and  infinite  impulse  response 
(HR)  types)  smear  the  boundries  between  (iie/a«(-  and  slow- 
phases  of  nystagmus  and  the  fixation  and  fast  components 
of  saccadic  EMs. 

We  have  adapted  a  nonlinear  smoothing  filter  (originally 
designed  to  optimise  edge  preservation  in  image  processing 
applications)  for  tile  amootiiinf  of  EM  signals.  This  filter  is 
called  a  Predictive  FIR-Median  Hybrid  (PFMH)  filter.  The 
PFMH  filter  operates  on  a  moving  window  of  data  samples 
centered  at  the  current  point  of  interest.  Several  predictive 
FIR  filters  are  applied  to  the  “upper”  and  “lower”  halves 
of  the  window  and  each  are  designed  to  predict  tiie  sample 
value  at  tiie  center  of  tile  windoir.  The  median  of  these  FIR 
filter  outputs  and  the  actual  center  data  sample  are  tsJcen 
as  the  PFMH  filter  output  for  each  window  position.  By 
properly  choosing  the  length  and  structure  of  the  FIR  sub- 
filters,  a  PFMH  filter  can  be  designed  to  smooth  a  bimodal 
EM  signal  without  blurring  the  boundries  between  the  two 
signal  components. 


Key  Words;  Non-Linear  Filters,  Order-Statistic  Filters, 
Eye  Movements,  Nystagmus,  Signal  Processing. 


John  D.  Enderle 

Division  of  Biomedical  Engineering 
Department  of  Electrical  and 
Electronics  Engineering 
North  Dakota  State  University 
Fargo,  ND  58105 

INTRODUCTION 

In  the  clinical  environment,  eye-movements  (EMs)  are  gen¬ 
erally  recorded  using  electrooculographic  (EOG)  methods. 
The  EOG  signal  results  from  the  nominal  0.4-0. 7  mv  poten¬ 
tial  difference  between  the  cornea  and  retina.  This  small  po¬ 
tential  (cornea  positive  with  respect  to  the  retina)  provides 
each  eyeball  with  an  electric  dipole  moment.  A  signal  pro¬ 
portional  to  the  angular  position  of  the  eyes  can  be  recorded 
by  placing  electrodes  at  the  outer  canthi  (outermost  edge)  of 
each  eye  and  amplifying  the  voltage  difference  between  the 
electrodes  [8].  The  resulting  EM  signal  has  a  scale  factor  of 
about  20pv/degree.  The  EM  signal  level  is  about  the  same 
magnitude  as  EMG  and  large  EEG  potentials.  Therefore, 
it  is  not  suprising  that  EM  signals  recorded  in  this  way  are 
often  corrupted  by  EMG  signals  from  the  facial  muscles  and 
sometimes  EEG  potentials.  It  has  been  standard  practice 
to  use  linear  low-pass  filters  to  remove  the  high-frequency 
EMG  noise  from  the  EM  recordings.  The  low-level,  low- 
frequency  EEG  noise  is  usually  tolerated. 

The  design  of  suitable  smoothing  filters  is  compUcated  by 
the  fact  that  some  EM  signals  are  bimodal.  Nystagmus 
consists  of  a  sequence  of  slow-  and  fast-phase  movements. 
The  slow-  and  fast-phases  of  nystagmus  are  generated  by 
different  neural  drcuits  and  have  different  dynamic  charac- 
terstics.  Ordinary  low-pass  filters  smear  the  boundries  be¬ 
tween  the  slow-  and  fast-phases,  and  therefore,  tend  to  dis¬ 
tort  the  signal.  Reducing  the  filter  cutoff  frequency  enough 
to  remove  the  noise  can  sometimes  result  in  severe  distor¬ 
tion.  The  same  applies  to  saccadic  EMs,  each  saccade  is 
separated  from  the  next  by  fixation  intervals.  The  signal 
characterics  during  fixation  are  different  than  during  a  sac- 
cade.  Linear,  low-pass  filters  blur  the  transition  between 
the  fixation  and  saccadic  signal  components. 


FILTER  DESIGN 

We  have  developed  a  method  for  smoothing  nystagmus  and 
saccadic  EM  signsJs  using  a  nonlinear  filter  called  a  Pre¬ 
dictive  FIR-Median  Hybrid  (PFMH)  filter  (4).  This  filter 
design  is  based  on  the  concept  that  both  nystagmus  and 
saccaxlic  EMs  can  be  modeled  accurately  as  a  sequence  of 
piecewise  continuous  polynomial  segments.  The  PFMH  fil¬ 
ter  is  designed  to  attenuate  signal  components  that  are  not 
consistent  with  this  signal  model.  Sharp  changes  between 
well  behaved  regions  (such  as  the  transition  between  slow- 
and  fast-phases  of  nystagmus)  are  consistant  with  the  signal 
model  and  are  preserved.  Oscillations  and  impulses  are  not 
consistent  with  the  signal  model  and  are  suppressed.  The 
PFMH  filter  smoothes  the  individual  signal  components  and 
preserves  the  boundries  between  them.  The  PFMH  filter 
acts  like  a  low-pass  filter  that  respects  the  natural  transi¬ 
tions  in  the  signal  and  preserves  them. 

Order  Statistic  Filters 

The  PFMH  filter  is  a  member  of  a  class  of  nonlinear  fillers 
called  Order-Statistic  (OS)  filters  [1,  6|.  OS  filters  operate 
on  a  moving  window  of  data  samples.  Generally,  the  win¬ 
dow  is  of  odd  length,  i.e.,  L  =  2N  -f  1.  The  data  samples  in 
the  window  are  rank  ordered  and  tl  j  rank-ordered  samples 
(order-statistics)  are  then  linearly  weighted.  The  linearly 
weighted  combination  of  order  statistics  constitute  the  fil¬ 
ter  output.  For  example,  the  standard  median  filter  is  an 
OS  filter  in  which  the  center  sample  of  the  ordered  window 
is  given  a  weight  of  I  and  the  other  samples  a  weight  of  0. 
Various  other  OS  filters  are  implemented  by  choosing  differ¬ 
ent  weighting  schemes.  The  OS  filter  concept  includes  the 
PFMH  filter,  a  design  developed  by  Heinonen  and  Neuvo 
(3,  4].  In  the  PFMH  filter,  the  data  samples  in  the  win¬ 
dow  are  first  processed  by  a  set  of  linear  FIR  predictive 
filters,  each  designed  to  predict  the  center  value  in  the  win¬ 
dow.  The  PFMH  filter  output  is  the  median  of  the  FIR 
filter  outputs  and  the  acutal  center  data  value.  The  dis- 
tintive  feature  of  all  OS  filters  is  the  rank-ordering  process, 
a  data  dependent,  nonlinear  operation.  OS  filters  are  non¬ 
linear  in  the  sense  that  superposition  does  not  apply.  If 
we  define  the  operation  of  OS  filtering  as  OS{  },  and  u(n) 
and  v(n)  are  different  time  varing  signals  then,  in  general, 
05{u(n)  -^  v(n)}  ^  (?5{u(n)}  -1-  05{v(n)}.  OS  filter*  'xe 
translation  and  offset  invariant  since,  05(4  •  tt(n)  -f-  H}  = 
A  •  OS{n(n)}  -t  B,  where  A  and  B  are  constants. 

Some  OS  filters,  such  as  the  median  and  PFMH  filters,  pos¬ 
ses*  a  set  of  signals  called  root  signals  that  are  invarient  to 
the  filter  [2,  4,  5,  7].  That  is,  toot  signals  will  pass  through 
the  filter  unmodified.  Also,  any  arbitrary  signal  will  be  re¬ 
duced  to  a  root  signal  after  a  finite  number  of  passe*  through 
the  filter. 


Smoothing  the  Signal 

One  would  expect  that  if  an  underlying  signal  were  a  root 
signal  for  a  given  OS  filter,  then  that  filter  could  extract 
the  signal  from  noise  without  distortion.  Repeated  filtering 
should  reduce  the  noise  contaminated  signal  to  a  root  signal, 
and  thereby  recover  the  signal  and  reject  the  noise.  Since 
the  number  of  toot  signals  for  a  pven  OS  filter  is  infinite, 
it  is  possible  that  the  filtering  process  might  not  converge 
to  the  desired  root  signal.  We  have  found,  in  practice,  that 
is  not  necessary  to  reduce  the  EM  signal  to  a  root  signal 
to  obtain  adequate  smoothing;  one  or  two  passes  through 
the  filter  provides  a  significant  noise  reduction.  In  any  case, 
repeated  filtering  generally  converges  to  a  root  signal  very 
close  to  the  desired  signal  in  about  five  passes,  or  so. 

Designing  the  PFMH  Filter 

The  PFMH  filter  operates  on  a  moving  window,  VV(n),  of 
data  samples,  x(y),  of  length  L  —  2N  -4-  !.  such  that 

W(t»)  =  {i!(n-  iV),...,x(n),...,*(n-i-  S)] 

The  data  samples  btfort  and  afitr  the  center  sample,  x(n), 
are  used  to  estimate  the  value  at  the  center  using  predictive 
FIR  substructures  of  the  form: 

Forward  predictor, 

=  E  •  x(n  -  i),  »  =  1, . . . ,  JV 

tsl 

Backward  predictor,  ' 

N 

*«(’»)  =  I2^rl»)  x(n  -t-t),  i  =  1,...,JV 


Where  if  (n)  and  iB(n)  aie  the  forward  and  backward  pre¬ 
dicted  values  and  fir(t)  is  the  coefficient  array  for  the  FIR 
predictor  of  order  r.  The  output  of  the  PFMH  filler  is  the 
median  of  the  predicted  values  and  the  data  value  x(n). 
Thus,  if  the  median  operation  ir  defined  a*  M£D{'},then 
the  filter  output  y(n)  is  given  by 

y(n)  =  M ED{i'f  {n), . . . ,  r(n), . . . ,  iB(n)} 

The  ee’-o-,  first-,  and  second-order  FIR  predictor  coeffidents 
as  given  by  Heinonen  and  Neuvo  (4)  are 


h^(i)  =  l/N,  i  =  1 . N 


h.(i)  = 


AN  -6i  + 2 
N{N  -  1) 


,  1,...,W 


.  ,  91V* -I- (9  -  36i)Ar -I- 30i^  -  18i  4- 6  . 


V 

. . ,  i » 


These  predictive  FIR  filter  subitrnctares  are  optimal  pre¬ 
dictors  for  rth  order  polynomials  corrupted  with  additive 
Gaussian  white  noise.  In  practice,  polynomials  up  to  about 
third-order  are  used.  We  have  been  using  four  predictive 
substructures  in  our  filters.  For  smoothing  nystagmus,  we 
use  three  first-order  predictors  of  different  lengths,  and  one 
second-order  predictor.  Saccadic  EMs  are  smoothed  using 
one  zero-order,  two  first-order,  and  one  second-order  predic- 
itor.  Several  predictors  of  different  length  and/or  order  are 
used  to  insure  having  at  least  one  that  works  well  for  each 
part  of  the  signed.  The  median  operation  selects  one  of  the 
predicted  values  (or  the  center  data  value)  as  the  filter  out¬ 
put. 

METHODS  AND  RESULTS 

We  recorded  EMs  using  standard  EOG  methods.  The  sig¬ 
nals  were  recorded  using  ailvei — silver  chloride  bitemporal 
electrodes.  The  amplified  electrode  potentials  were  then 
passed  through  a  30-Hz,  low-pass,  4-pole  Butterworth  filter 
SLnd  digitized  at  a  sampling  rate  of  128  Hz.  A  PFMH  fil¬ 
ter  was  designed  for  nystagmus  smoothing  using  three  first- 
order  predictors  of  sizes,  IV  =  6,  iV  =  15  and  N  =  24,  and  a 
second-order  predictor  of  size,  N  =  24. 

Fig.  1  shows  the  results  of  successively  filtering  a  noisy  nys¬ 
tagmus  signal  five  times  using  the  PFMH  filter.  After  the 
second  pass,  the  noise  is  nearly  gone.  After  five  passes,  the 
signal  is  very  near  a  root  signal  and  almost  noise  free.  The 
transitions  between  the  slow-  and  fast-phases  are  actually 
enhanced  by  successive  passes  of  the  PFMH  filter.  This 
sharpening  cf  the  “cornets”  on  the  waveform  is  a  result  of 
the  root  structure  of  the  filter.  Repeated  filtering  using  a 
I'near  low-pass  filter  would  increasingly  blur  the  tnuisitions. 

In  order  to  provide  a  comparison  of  the  PFMH  filter  to  a 
standard  low-pass  filter,  we  corrupted  an  idealized  nystag¬ 
mus  signal  by  adding  30-Hz  low-pass  filtered  noise.  The 
“ideal”  signal  consisted  of  a  sequence  of  slow-phase  seg¬ 
ments  represented  by  second-order  polynomial  segments, 
separated  by  linear  high-velodty,  fast-phase  segments.  We 
then  processed  the  corrupted  signal  with  each  filter  sepa¬ 
rately.  The  signal  was  passed  through  each  filter  once.  The 
filtered  signals  were  then  compared  to  the  uncorrupted  origi¬ 
nal.  The  results  are  illustrated  in  Fig.  2.  The  mean-squared 
difference  (MSD)  for  each  filter  was  computed  This  was 
done  by  calculating  the  sum  of  the  squared  differences  be¬ 
tween  each  filtered  signal  and  the  uncorrupted  signal.  The 
MSD  for  the  PFMH  filter  was  6755  and  the  MSD  for  the 
low-past  filter  was  8574.  Thus,  the  linear  low-pats  filter  had 
a  27%  larger  MSD  than  the  PFMH  filte-  P.»e»v'*!ntio"  of 
Fig  ?  discloses  that  one  pass  of  the  PFMH  filter  effectively 
smooths  the  slow-phase  components  and  yet  preserves  the 
slow-  to  fast-phase  transitions.  By  contrast,  the  low-pass  fil¬ 
ter  allows  considerable  low-frequency  noise  to  pass  through 
while  blurring  the  signal  transitions. 


DISCUSSION 

In  most  applications,  the  standard  linear  low-pass  filter  is 
an  efficient,  effective  smoothing  filter.  This  is  particularly 
true  for  signals  generated  by  stationary  processes  when  the 
spectrum  of  the  corrupting  noise  does  not  significantly  over¬ 
lap  the  signal  spectrum  and  the  noise  has  a  Gaussian  am¬ 
plitude  distribution.  Many  practical  data  smoothing  appli¬ 
cations  meet  these  conditions — at  least  approximately.  In 
other  instances,  the  signals  are  clearly  multi-modal  and  are 
not  generated  by  stationary  processes;  this  occurs  for  many 
biological  signals.  These  nonstationary  signals  can  often  be 
smoothed  more  effectively  using  a  nonlinear  filter  such  as 
the  PFMH  filter. 

For  the  PFMH  filter  to  be  effective,  the  basic  underlying 
signal  to  be  smoothed  must  be  a  root  signal  for  the  filter. 
This  requirement  imposes  a  limitation  for  the  method  since 
only  a  limited  number  of  PFMH  filter  designs  have  been  dis¬ 
covered.  To  date,  PFMH  filters  have  been  designed  using 
linear  FIR  substructures  and  having  root  signals  that  are 
piecewise  continuous,  low-order  polynomials. 

Since  the  PFMH  filter  is  nonlinear  it  does  not  have  a  trans¬ 
fer  function  or  any  other  frequency  domain  description.  The 
principle  of  superposition  does  not  apply,  so  knowing  the 
sine-wave  response  does  not  provide  any  information  as  to 
how  complex  signals  will  be  processed.  Therefore,  there  is 
no  simple  way  to  describe  the  “bandwidth”  or  “frequency 
response”  of  the  PFMH  filter. 

Selecting  the  best  filter  parameters  (i.e.,  length,  order,  and 
number  of  FIR  substructures)  is  difficult.  There  are  no  de¬ 
sign  formulas  to  provide  optimum  parameters  for  any  given 
application.  The  best  epproach  is  to  use  a  Monte  Carlo  pro¬ 
cedure  to  search  a  range  of  parameters  using  a  large,  rep¬ 
resentative  data  set.  The  data  set  may  consist  of  simulated 
data.  The  parameter  values  that  perform  best  according 
to  some  criteria,  such  as  MSD,  will  generally  work  well  for 
actual  data. 

Biological  signals  are  often  corrupted  by  non-Gaussian  and 
sometimes  non-stationary  noise  sources.  For  example,  the 
typical  electrocardiogram  (ECG)  signal  has  a  non-Gaussian 
amplitude  distribution,  mainly  due  to  the  large  QRS  com¬ 
plex.  In  cases  where  the  ECG  signal  is  a  noise  source  (as  in 
EM  recording),  the  corrupting  noise  is  clearly  non-Gaussian. 
Since  the  currently  available  t»FMII  filters  ha'^e  HTi»‘aT  sub¬ 
structures,  they  are  optimum  smoothing  filters  for  signals 
corrupted  by  additive  Gaussian  noise.  These  filter  designs 
are  less  effective  for  removal  -'f  nc%  oatusian  ntnse  such  as 
Impulse  uoise  or  noise  having  skewed  or  heavy-tailed  am¬ 
plitude  distributions.  Other,  perhaps  nonlinear,  predictive 
FIR  substructures  could  possibly  be  used  to  generate  new 
PFMH  filters  having  new  root  signal  structures  and  im¬ 
proved  performance  for  non-Gaussian  noise.  We  are  plan¬ 
ning  to  explore  several  filter  designs  of  this  type. 
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FIGURE  1.  The  PFMH  Alter  is  used  to  Alter  e  perticuierly  noisy  segment  of  nysteg- 
mus.  Top  trace  is  the  noisy  input  signal,  the  Ave  following  traces  are  successively 
Altered  versions  of  the  signal.  Each  pass  through  the  Alter  brings  the  nystag¬ 
mus  signal  closer  to  a  root  signal.  One  or  two  passes  produce  a  signiAcant  noise 
reduction. 
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FIGURE  2.  An  nrtiflcinl  nystagmus'like  signal  was  created  using  second*  order  slow- 
phase  components  and  flrst*order  fast-phase  components.  This  signal  is  a  root 
signal  for  the  PFMH  Alter.  The  ideal  signal  (top  trace)  was  corrupted  by  adding 
Gaussian  noise  and  the  corrupted  signal  (second  trace)  was  then  Altered  separately 
by  the  PFMH  Alter  and  a  linear,  10  Hs  low-pass  Alter.  An  error  signal  was 
calculated  by  subtracting  the  Altered  signals  from  the  uncormpted  signal.  The 
PFMH  Altered  signal  and  the  corresponding  error  signal  (multiplied  by  6)  are 
shown  in  the  center  two  traces.  The  iast  two  traces  show  Altered  signal  and  the 
X5  error  signal  for  the  10-Hs  low-pass  Alter. 


